samp_n <- 100
rand_samp_data <- read.csv(file.path(dir_M, "git-annex/impact_acceleration/stressors/sst/int/rand_sample_diff.csv"))
for(i in 1:samp_n){ # i = 2
plot_sst <- ggplot(dplyr::filter(rand_samp_data, sample_id==i), aes(y = sst_anom, x = year)) +
geom_point(size = 2) +
geom_line() +
stat_smooth(method=lm, se=FALSE, color="red")
plot(plot_sst)
}